Hidden magnetic frustration by quantum relaxation in anisotropic Nd-langasite 



00 

o 
o 

(N 
C 

3 



V. Simonet, 1 R. Ballou, 1 J. Robert, 1 B. Canals, 1 F. Hippert, 2 P. 
Bordet, 1 P. Lejay, 1 P. Fouquet, 3 J. Ollivier, 3 and D. Braithwaite 4 

'institut Neel, CNRS/UJF, B.P. 166, 38 042 Grenoble Cedex 9, France 
2 LMGP (CNRS), Grenoble Institute of Technology, 38016 Grenoble, France 
Institut Laue-Langevin, BP 156, 38042 Grenoble Cedex 9, France 
4 CEA-Grenoble, INAC/SPSMS/IMAPEC, 17 rue des Martyrs, 38054 Grenoble cedex, France 

(Dated: June 13, 2008) 

The static and dynamic magnetic properties of the Nd3GasSiOi4 compound, which appears as the 
first materialization of a rare-earth kagome-type lattice, were re-examined, owing to contradictory 
results in the previous studies. Neutron scattering, magnetization and specific heat measurements 
were performed and analyzed, in particular by fully taking account of the crystal electric field 
effects on the Nd 3+ ions. One of the novel findings is that the peculiar temperature independent 
spin dynamics observed below 10 K expresses single-ion quantum processes. This would short-circuit 
the frustration induced cooperative dynamics, which would emerge only at very low temperature. 
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The kagome antiferromagnet is the prototype of geo- 
metrically frustrated systems in 2 dimensions. It is pre- 
dicted that its ground state, in the classical case, remains 
disordered and is highly degenerate The investiga- 
tion of such a spin-liquid state is still impeded by the 
few available experimental materializations. Various sec- 
ondary perturbations actually lead, in most known com- 
pounds, to long-range order (LRO) or spin glass phase. 
This occurs, though, at a temperature well below the 
onset of magnetic correlations and is often accompanied 
by a dynamical behavior persisting below the transition 
temperature. These dynamics are characterized by a 
slowing down of the magnetic fluctuations which rounds 
off onto a relaxation plateau Q . Although rather generic, 
the origin of this peculiar relaxation is still debated. 

In all the kagome compounds so far studied the mag- 
netic entities that form the frustrated lattice of corner- 
sharing triangles are 3d transition metal ions with weak 
magneto-crystalline anisotropy. Inversely, in the 3 di- 
mensional compounds of the pyrochlore family, 4f rare- 
earth ions with strong magneto-crystalline anisotropy oc- 
cupy a frustrated lattice of corner-sharing tetrahedra. 
This leads to a wider variety of magnetic behaviors. Most 
striking, for instance, was the discovery of the spin ice 
state in Ho2Ti207, where the frustration arises from 
the combined effects of multi-axial anisotropy and fer- 
romagnetic interactions 0]. It was expected, in con- 
trast, that the multi-axial anisotropy would release the 
geometric frustration of antiferromagnetic interactions in 
Tb2Ti207, but this compound exhibits no LRO and, 
instead, shows intriguing spin dynamics well below the 
Curie- Weiss (CW) temperature [5j, |6j . It was recently ar- 
gued that this compound would form a quantum spin ice 
from dynamical frustration induced through crystal field 
(CF) excitations and quantum many-body effects Q. 

The lately discovered NdsGasSiOu (NGS) compound, 
which belongs to the langasite family, provides with the 



unique opportunity to study the role of a strong magneto- 
crystalline anisotropy in a geometrically frustrated 2 di- 
mensional lattice. In this compound the rare- earth Nd 3+ 
ions fully occupy a lattice showing the kagome corner- 
sharing triangles connectivity if only first neighbor inter- 
actions are taken into consideration Q . 

The magnetic properties of NGS powder samples and 
single crystal were investigated in the 1.6—400 K temper- 
ature range 8]. Neutron diffraction and magnetization 
measurements on the powders confirmed that NGS does 
not order in a Neel phase nor shows any signature of a 
spin glass behaviour. Magnetization measurements on 
the single-crystal (cf. Fig. HJ revealed a strong magneto- 
crystalline anisotropy associated with CF effects on the 
ground multiplet J=9/2 of the Nd 3+ ions. A mean field 
quantitative analysis of the magnetic susceptibility at 
high temperature taking into account only the quadrupo- 
lar contribution to the CF yielded a CW temperature 
6 = —52 K, leading to anticipate rather strong antiferro- 
magnetic interactions between the Nd moments. These 
were described as coplanar rotators in the kagome planes. 
A change in the anisotropy occurs around 33 K, below 
which the c-axis (perpendicular to the kagome planes) 
becomes the macroscopic magnetization axis. 

Recently, a neutron scattering investigation of NGS 
at very low temperature was reported [9(. A weak ring 
of diffuse magnetic scattering was observed at 46 mK, 
whose momentum-transfer (Q) dependence was found 
in agreement with calculated magnetic scattering maps 
from spin- liquid models dfj]. It was inferred from addi- 
tional neutron spin-echo (NSE) experiments, that the as- 
sociated magnetic state would be highly dynamical. The 
authors also claimed to have evidenced a magnetic field 
induced phase transition featured by a plateau of reduced 
magnetization, approximately half the moment value of 
the free Nd 3+ ion, in agreement with magnetization mea- 
surements (1.6 [Ib at 1.6 K) Q. Another recent investiga- 
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FIG. 1: (Color online) Time-of-flight spectra, Q-integrated in 
the range [0.34, 2.55] A -1 , a) Quasielastic magnetic scatter- 
ing fitted from 75 to 300 K by a Lorentz function multiplied 
by the detailed balance factor. The fluctuation time obtained 
from the inverse HWHM is shown in the inset together with 
an Arrhenius fit (r = r exp(AE/k B T) with AE=133 K and 
to=6.7 10~ 13 s). b) Two non-dispersive crystal-field level 
transitions pointed out by the arrows. Note that the feature at 
~ 23 meV was attributed to phonons from the Q-dependence. 



tion by //SR and NMR spectroscopy [ljj established that 
the magnetic fluctuations slow down as the temperature 
is decreased to reach a relaxation plateau below 10 K, 
interpreted as a signature of magnetic frustration in this 
cooperative system. Some usual signatures of frustration 
therefore appeared to have been observed in NGS, mag- 
netic spatial correlations and a fluctuating ground state, 
but inconsistencies soon emerge from close inspection. In 
reference 0, no magnetic signal is reported at 46 mK out- 
side an energy window corresponding to fluctuating times 
larger than 5 1CP 11 s but NSE experiments at the same 
temperature appears to suggest a relaxation faster than 
4 10~ 12 s. In addition, this fast relaxation suggested 
by the NSE experiment is incompatible with the char- 
acteristic fluctuation time of the /iSR-probed relaxation 
plateau, r « 2.5 10~ 10 s. This motivated our exper- 
imental re-investigation of the spin dynamics in NGS, 
presented below, as well as a full characterization of the 
single-ion anisotropy through a complete CF analysis. 

Powders of Nd-langasite, synthesized as in ref. d, were 
used for magnetization and specific heat measurements 
and neutron scattering experiments. The spin dynamics 
were investigated via time-of-flight (ToF) and NSE spec- 
troscopy. The ToF experiments were performed at ILL 
on the IN4 spectrometer at a wavelength of 1.1 A and 
on the IN5 spectrometer with wavelengths between 8.5 
and 2.25 A. The ToF spectra reveal a quasielastic (QE) 
signal, displayed in Fig. [TJ whose magnetic origin was 
inferred from its Q-dependence (slight decrease follow- 
ing the Nd 3+ squared form factor). The half width at 
half maximum (HWHM) of the Q-integrated spectra de- 
creases from 300 K to 75 K. An Arrhenius law can be 
forced to account for the temperature dependence of the 
relaxation time r = S/HWHM. This yields an energy 
barrier of the order of 130 K for a thermally activated 



process. At 2 K, no QE signal can be detected any more 
within the instrumental resolution energy window. 

In order to extend the time window for the obser- 
vations of the dynamics, a NSE experiment was per- 
formed on the IN11C spectrometer at ILL using a wave- 
length of 5.5 A. XYZ polarization analysis was performed 
to experimentally select the magnetic part of the sig- 
nal. The instrumental resolution was corrected with the 
purely elastic magnetic signal of a reference compound. 
The accessible time window of the NSE technique over- 
laps the ToF range for short times and the fj.SK range 
for long times. The normalized intermediate functions 
S(Q, t) I S(Q, 0), measured between 1.4 and 100 K, can be 
fitted by stretched exponential line-shapes exp(— (t/r)^) 
(cf. Fig. [2|). A f3 exponent ^ 1 denotes a non- unique 
relaxation time. It increases gradually from a constant 
value of 0.5 between 1.4 and 10 K to 0.8 at 40 K, indicat- 
ing a change in the dynamics. The small temporal range 
with non-zero signal at higher temperatures prevented 
from obtaining a reliable value of (5 which was therefore 
fixed equal to 0.8 between 60 and 100 K. Above 40 K, 
the thermal variation of the relaxation time agrees with a 
thermally activated process (ln(r) oc 1/T). We get an en- 
ergy barrier very close to the one extracted from the ToF 
QE signal (cf. Fig. [2|) [13]. Below 10 K, the relaxation 
time does no more vary with temperature, revealing a re- 
laxation plateau with a characteristic time r « 3 10~ 10 s, 
in good agreement with the /iSR results [ll[ (and at vari- 
ance with the NSE data of ref. l9|). An important finding 
of our NSE experiment is that the magnetic relaxation 
does not show any detectable Q-dependence in the range 
[0.3 A" 1 , 1.7 A -1 ], pointing out the absence of spatial 
correlation down to 1.4 K since the frequency width is 
inversely proportional to the Q-dependent susceptibility 
in correlated paramagnet (l3| . The spatial correlations 
observed at 46 mK Q thus would occur in a temperature 
range much below the onset of the relaxation plateau. 

Additional insights about the magnetism of NGS were 
gained from the ToF scattering at higher energy transfer, 
specific heat measurements and investigations of magnet- 
ically diluted samples. Two CF level transitions can be 
identified, through their Q dependence in the ToF spec- 
tra, around 95 K (also reported in ref. d) and 380 K 
(cf. Fig. [I]). They however are very broad in energy, with 
a width at least one order of magnitude larger than the 
instrumental resolution. Quantitative estimations of CF 
effects suggest that this may be ascribed to the distri- 
bution of the charge environments associated with the 
Si 4+ /Ga 3+ mixed site which also is responsible for 
the NMR line broadening reported in ref.lll|. The specific 
heat measurements were performed down to 400 mK in a 
Quantum Design PPMS. The magnetic contribution was 
determined after subtraction of the specific heat mea- 
sured on the isostructural non-magnetic LaaGasSiOu 
compound. It exhibits a broad bump ranging from 
«100 K to above 300 K. A strong rise is noticeable below 
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FIG. 2: (Color online) Neutron spin-echo normalized interme- 
diate functions fitted with stretched exponentials (solid lines) 
for several temperatures (a), and for several Q ranges at 1.4 
K (b). The variation of the relaxation time over the inverse 
temperature is displayed with an Arrhenius fit (solid line) at 
high temperatures (AE=135 K, to — 1.2 1CF 12 s) and a dashed 
line underlying the relaxation plateau below 10 K (c). 

2 K as the temperature is decreased down to the lowest 
values (cf. Fig. [3]). We shall discuss this below. Also 
meaningful were the magnetization measurements on 
the magnetically diluted (Ndo.oiLao.99)3Ga5SiOi4 com- 
pound, which were performed in order to isolate the in- 
fluence of the magneto-crystalline anisotropy. We found 
out, to our surprise, that its magnetic susceptibility does 
experimentally not differ from that of NGS (cf. Fig. [3} , 
demonstrating that the exchange interaction between the 
Nd + moments in NGS is certainly much smaller than 
initially estimated [8|, at least not greater than a few K. 

The above results were confronted to calculations of 
the CF effects. In NGS, there are 3 distinct Nd 3+ ions 
per unit cell, each lying on a different 2-fold axis, inter- 
related through the 3-fold c axis of the crystal [8[ , thus 
giving rise to multi- axial low symmetry CF potentials. 
The associated Hamiltonian for each Nd 3+ ion in the 
ground multiplet J=9/2 can be expressed as H CF — 

EU^-P A2 2 q P ° 2 4 with K q = (^)*. in terms of 

O™ Racah operators and with respect to a local frame 
where the 2-fold axis is the quantization axis for the angu- 
lar momentum Ignoring exchange interactions this 
provides with an energy spectrum which can be compared 




100 200 300 1 2 3 



T (K) T(K) 

FIG. 3: (Color online) a) Inverse susceptibility of NGS and 
of a magnetically diluted sample, b) NGS low temperature 

exp ( k AE r ) 

specific heat fitted by R(]^p|0 2 7j[^; ( — B At; p f° r a Schottky 

anomaly produced by a two-level system with R the gas con- 
stant and AE=235 mK. 



with the transitions detected in the ToF spectra and can 
be used to calculate the specific heat. The magnetiza- 
tion and susceptibility are obtained from the eigenfunc- 
tions and must be averaged over the 3 Nd 3+ ions. This 
forced us to work in a common global frame with the c 
axis as the quantization axis and to perform appropriate 
rotations of the Racah operators. It was illusory to try 
to determine the 9 independent A 2 ^ p coefficients at once. 
Thus, in a first step the A\ p jA^ v ratios were calculated in 
the point charge model from the first coordination shell 
of eight 2 ~ (cf. Fig. [4]). Qualitatively good results were 
already obtained, reproducing in particular the crossing 
of the susceptibility for applied fields parallel and perpen- 
dicular to c. This crossing is due to the 4 th order terms of 
H CF . The A 2 2 q p coefficients 

were, in a second step, slightly 
varied in order to better account for all the experimental 
constrains (details of the calculations will be given else- 
where) . With the best set of A^ p parameters, the energy 
spectrum consists of 5 degenerate doublets, the first two 
being separated from the ground state doublet by 90 K 
and 324 K respectively, in rather good agreement with 
the ToF spectra considering the width of the measured 
excitations. The ground state in the global frame with 
the c quantization axis is a combination of | ± 7/2 >, 
| ±5/2 > and | ±3/2 > states with a maximum weight 
for the | ±5/2 > state. The good agreement between the 
calculated and measured magnetic isotherms and suscep- 
tibility is shown in Fig. |H A first observation is that the 
low temperature reduced moment must be attributed to 
a CF effect, the so-called field-induced magnetic phase 
0, [ll[ being hence a polarized paramagnet, at least at 
1.6 K, as in Tb2Ti 2 07 A second observation is that 
the susceptibility reaches the high temperature behaviour 
above 2000 K. This regime occurs thus at a temperature 
much higher than the one associated to the highest crys- 
tal field level energy (665 K). This is a consequence of 
the large 6 th order terms, shifting the inverse susceptibil- 
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FIG. 4: (Color online) Full CF calculation (solid lines) of the 
inverse susceptibility (main frame) and magnetic isotherms 
(inset) with the magnetic field parallel (red) or perpendicu- 
lar (blue) to the kagome planes compared to single-crystal 
measurements (circles) Q|. The dashed lines are high tem- 
perature expansion calculations. The three Nd 3+ are shown 
in their local frame with their first coordination O 2- shells. 



ities upward, which explains why the previous estimation 
of the exchange interactions went erroneous despite the 
almost linear variations of the inverse susceptibility [|| . 

The CF analysis confirms the smallness of the exchange 
interaction in consistency with the absence of spatial cor- 
relations down to 1.4 K. This suggests a single-ion origin 
for the measured magnetic relaxation both in the vary- 
ing and independent temperature regimes, as probed by 
the neutron spectroscopy. The first excited doublet of 
the calculated CF spectrum lies at an energy in qualita- 
tive agreement with the barrier extracted from the ToF 
and NSE measurements at temperatures above 40 K. The 
thermally activated dynamics in this temperature range 
are most probably mediated by the spin-phonon interac- 
tion [l7j. On decreasing the temperature, this progres- 
sively slows down, but below 10 K is short-circuited by 
another process with a temperature independent relax- 
ation time. This strongly suggests quantum relaxation 
phenomena 18, IT^ | and calls to mind the dynamical 
cross-over observed in spin ice compounds [20| . The ori- 
gin of this relaxation is still unclear but could be related 
to the only feature not accounted for by the CF calcula- 
tion: the rise of the specific heat at very low temperature. 
A splitting of the ground state doublet was suggested in 
ref. © from extrapolation of Tof data to zero-field under 
magnetic field. Assuming this suggestion, we fitted the 
specific heat by a T~ 2 law in agreement with the tail of a 
Schottky anomaly associated with a two-level system (cf. 
Fig. [3]). This gives an energy gap of «235 mK, in good 
agreement with the reported zero-field extrapolated split- 
ting Although the cause of this splitting is uniden- 
tified, we know that it should necessarily break the time 



reversal symmetry, since the Nd 3+ are Kramers ions. An 
interaction involving odd powers of the ladder operators 
of the Nd 3+ angular momentum would do, by providing 
with non-diagonal matrix elements mixing the two or- 
thogonal states of the doublet. The consequences of this 
low temperature dynamics on the cooperative magnetic 
phase, expected when the exchange interactions and the 
geometric frustration become significant, can now be ad- 
dressed. A cooperative dynamical behaviour presenting 
longer fluctuation times than the fluctuations created by 
the proposed single-ion quantum process would indeed 
be hindered. This, thus, queries the nature of the phase 
associated to the ring-like Q-distribution at 46 mK . 

In conclusion, the experimental re- investigations of the 
kagome-like rare-earth NGS compound combined with 
CF analysis have allowed featuring its magnetism on firm 
grounds, solving previous inconsistencies. The interpre- 
tation of the relaxation plateau of the moment dynamics 
in terms of cooperative processes is ruled out. Its origin 
must rather be sought from single-ion quantum processes. 
The effects of the geometrical frustration then would be 
hidden emerging only at very low temperature. 
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